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Principal  Component  Analysis  (PCA)  has  been  used  to  process  the  data  obtained  from  the  batch  pyrolysis 
of  high  density  polyethylene  (HDPE)  in  order  to  recognize  trends  in  the  formation  rates  of  the  large 
number  of  components  in  the  product  stream.  Two  different  criteria  for  component  grouping  have  been 
compared,  i.e.,  according  to  molecular  size,  and  the  nature  of  C— C  bond.  The  choice  of  the  suitable  crite¬ 
rion  is  crucial  in  order  to  establish  the  kinetic  scheme,  and  PCA  has  proven  to  be  a  very  helpful  tool.  The 
multivariate  analysis  applied  to  the  lumps  obtained  according  to  the  molecular  size  of  the  components 
provides  very  useful  information,  which  allows  proposing  similar  kinetic  schemes  to  those  in  the  litera¬ 
ture,  but  with  significant  details  involved  in  the  definition  of  the  lumps  and  their  distribution  between 
primary  and  secondary  reactions.  Furthermore,  although  there  are  substantial  differences  in  the  forma¬ 
tion  rate  of  single-ring  aromatics  (which  are  particularly  affected  by  temperature),  PCA  carried  out  by 
grouping  the  products  according  to  the  nature  of  their  C— C  bond  does  not  provide  relevant  information 
concerning  differences  in  the  performance  of  the  main  products  (aliphatic  hydrocarbons). 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

According  to  2013  Annual  Report  of  the  European  Association  of 
Plastics  Manufacturers  [1]  plastic  production  in  2012  amounted  to 
57  M  tonnes  in  Europe  and  288  M  tonnes  worldwide.  Thanks  to  a 
growing  public  awareness  and  the  consequent  improvements  in 
plastic  waste  management  policies,  landfill  has  decreased  in 
EU-27  over  the  last  5  years.  In  fact,  the  recovered  fraction  is  now 
higher  than  the  amount  sent  to  disposal/landfill  sites  since  2008. 
Recycling  (both  mechanical  and  chemical)  is  the  preferred  option 
for  waste  plastic,  as  the  product  shelf  life  is  increased.  Neverthe¬ 
less,  less  than  25%  of  the  25.1  M  tonnes  of  waste  plastic  generated 
in  Europe  in  2011  was  suitable  for  recycling,  which  is  due  to  the 
production  of  composites  and  specially  designed  plastics,  the  use 
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of  new  methods  that  improve  the  mechanical  and  superficial  prop¬ 
erties  of  the  plastics,  and  the  high  cost  of  separation  processes.  The 
rest  of  the  recovered  fraction  is  intended  for  energy  recovery 
(14.9  M  tonnes  recovered,  i.e.,  59%  of  waste  plastic). 

The  polyolefins  (HDPE,  LDPE,  LLDPE,  PP)  and  the  polyaromatics 
(PS,  EPS)  account  for  more  than  55%  of  the  European  plastic  pro¬ 
duction.  Although  they  have  similar  physical  and  chemical  proper¬ 
ties  and  their  individual  separation  is  complicated  and  cost 
intensive,  this  fraction  could  be  easily  isolated  from  other  MSW / 
plastic  fractions  using  commercially  available  density-based  sepa¬ 
ration  methods  [2].  Although  this  fraction  would  not  be  valid  for 
mechanical  recycling  owing  to  its  composition,  properties,  and 
thermal  behavior,  it  can  be  upgraded  by  chemical  recycling. 
Thermochemical  methods  would  allow  the  (re)generation  of 
hydrocarbon  feedstock,  and  therefore  mitigate  crude  demand  and 
fill  in  the  gap  between  waste-to-energy  and  mechanical  recycling 
[3,4  .  Pyrolysis  is  an  attractive  option  to  obtain  hydrocarbons  from 
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plastic  wastes,  as  has  been  reported  in  the  literature.  Aromatic 
yields  above  60  wt%  have  been  obtained  in  the  pyrolysis  of  polysty¬ 
rene  [5,6]  and  rejected  streams  of  packaging  plastic  wastes  from 
sorting  plants  [7];  more  than  55  wt%  of  light  olefins  in  the  catalytic 
pyrolysis  of  HDPE  [8-10];  and  yields  of  gasoline  above  55  wt%  in 
the  pyrolysis  of  HDPE  with  a  spent  FCC  catalyst  [11].  HDPE  has 
been  chosen  in  many  of  these  studies  as  a  plastic  waste  sample, 
because  the  behavior  of  the  rest  of  the  polyolefins  (LDPE,  LLDPE, 
ULDPE,  PP)  is  similar  [12-14]. 

These  results  illustrate  the  potential  pyrolysis  has  to  become  the 
plastics  chemical  recycling  (tertiary)  technology  par  excellence  and 
allow  for  optimism  in  the  search  for  solutions  to  the  pressing  prob¬ 
lem  involving  the  management  of  this  waste.  However,  scaling-up 
the  pyrolysis  reactor  from  lab  to  industrial  size  is  not  easy,  since 
there  is  an  insufficient  understanding  of  the  underlying  reaction 
pathways,  thereby  preventing  a  quantitative  prediction  of  the  full 
product  distribution  [15].  In  fact,  although  there  are  many  kinetic 
studies  regarding  the  pyrolysis  and  catalytic  pyrolysis  of  polyolefins, 
most  of  them  have  been  carried  out  in  thermogravimetric  analyzers 
[10,13,16-21  with  minor  samples  and  small  particles;  or  with 
devices  that  do  not  operate  under  conditions  equivalent  to  those 
of  industrial  units  [22,23];  that  is,  they  do  not  operate  under  ultra¬ 
fast  heating  and  continuous  plastic  waste  feed  and  product  stream. 
To  sum  up,  the  results  of  these  studies  should  be  taken  with  caution 
for  their  application  in  industrial  processes  [18  .  Kinetic  knowledge 
is  essential  for  the  determination  of  the  optimum  conditions  for 
maximizing  a  given  product  [24]  and  for  the  design  and  simulation 
of  commercial  reactors  [25].  Thus,  kinetic  studies  are  especially  rel¬ 
evant  when  they  have  been  carried  out  in  reactors  designed  to  over¬ 
come  the  current  limitations  for  scaling  up  industrial  ovens  [26-28]. 
These  reactors  should  have  the  same  configuration  as  commercial 
designs,  and  product  formation  should  be  modeled  according  to 
kinetic  schemes  that  adequately  predict  the  result  of  the  complex 
reactions  occurring  in  the  thermal  degradation  of  the  large  molecu¬ 
lar  chains  of  the  polymers. 

Regarding  the  kinetic  schemes  for  polyolefin  pyrolysis, 
Westerhout  et  al.  [29]  proposed  a  pioneering  model  in  which  the 
primary  devolatilization  reaction  yields  mainly  intermediate 
wax-like  products  that  are  cracked  to  smaller  hydrocarbons  fol¬ 
lowing  secondary  reactions.  These  authors  stated  that  low  molec¬ 
ular  weight  alkanes  and  alkenes  are  thermodynamically  unstable 
at  high  temperatures,  being  converted  into  aromatic  compounds 
(tertiary  reactions)  and,  in  the  case  of  sufficiently  high  residence 
times,  coke  could  be  formed.  In  a  subsequent  work,  Elordi  et  al. 
[30]  observed  that  this  model  does  not  effectively  predict  the 
results  obtained  in  a  conical  spouted  bed  reactor  in  the  450- 
715  °C  range.  This  disparity  between  the  aforementioned  models 
could  be  related  to  the  lower  amount  of  secondary  and  tertiary 
reactions  that  take  place  in  this  type  of  reactor  due  to  the  very 
short  residence  times  attained.  These  authors  tested  up  to  eight 
new  kinetic  schemes.  The  best  fit  was  achieved  with  a  model  in 
which  polyethylene  degrades  to  give  gas,  liquid,  and  wax;  with 
the  latter  undergoing  secondary  reactions  to  give  liquid  and  aro¬ 
matics,  which  in  turn  and  at  high  temperatures,  produces  more 
char.  In  a  more  recent  work,  Huang  et  al.  [28  proposed  a  mecha¬ 
nistic  scheme  including  the  main  features  and  kinetic  reaction 
schemes  for  polymer  degradation  over  cracking  catalysts,  and  dis¬ 
cussed  the  role  of  the  carbenium  ion  in  catalytic  cracking  chemis¬ 
try.  This  proposal  does  not  group  the  components  as  determined 
by  their  molecular  size  or  state  at  ambient  temperature  (gas, 
liquid,  wax),  but  puts  forward  another  criterion  based  on  the  nat¬ 
ure  of  the  chemical  bonds  that  are  characteristic  of  the  different 
molecules.  Hence,  products  are  lumped  into  paraffins,  olefins, 
and  aromatics. 

Selecting  the  grouping  criterion  that  best  represents  the  multi¬ 
ple  cracking  reactions  the  polymer  undergoes  is  a  key  factor  for  the 


proposal  and  validation  of  a  kinetic  scheme.  Accordingly,  analysis 
tools  are  required  to  objectively  and  quantitatively  elucidate  the 
benefits  of  one  criterion  over  another.  In  this  paper,  Principle 
Component  Analysis  (PCA)  is  proposed  in  order  to  select  the  best 
criterion  for  lumping  the  reaction  products,  given  that  this  mathe¬ 
matic  technique  of  multivariate  analysis  is  suitable  for  recognizing 
patterns  in  complex  problems  [31-34  ,  and  its  potential  has  been 
proven  in  the  kinetic  study  of  waste  material  pyrolysis  [35  . 

2.  Experimental 

Pyrolysis  experiments  have  been  carried  out  in  the  550-650  °C 
range  in  a  conical  spouted  bed  reactor  described  elsewhere  [36,37]. 
The  raw  material  is  HDPE  manufactured  by  Dow  Chemical 
(Tarragona,  Spain),  with  its  main  properties  being  as  follows:  aver¬ 
age  particle  size:  4  mm;  molecular  weight  (Mw):  46,200  g/ mol; 
polydispersity:  2.89;  density:  940  kg/m3;  and  higher  heating 
value:  46,230  kj/kg. 

Runs  have  been  carried  out  in  batch  mode.  The  reactor  bed  con¬ 
sists  of  30  g  of  sand  (particle  size  between  0.63  and  1  mm).  In  all 
cases  the  fluidizing  agent  was  N2,  with  a  velocity  of  1.2  times  of 
that  corresponding  to  the  minimum  spouting  regime.  1  g  of  HDPE 
is  fed  into  the  reactor,  and  during  each  run,  gas  samples  have  been 
taken  from  the  product  stream  at  given  reaction  times,  in  order  to 
carry  out  on-line  monitoring  by  gas  chromatography.  This  analysis 
method  allows  monitoring  the  evolution  over  time  of  all  the  vola¬ 
tile  products  with  a  molecular  size  smaller  than  20  carbon  atoms. 
The  remainder  of  products  correspond  to  heavy  fractions,  including 
light  waxes,  with  boiling  points  between  343  °C  and  525  °C,  and 
heavy  waxes,  with  higher  molecular  weights  and  boiling  tempera¬ 
tures  [38].  In  view  of  their  size,  the  characterization  of  waxes 
requires  more  specific  analytical  techniques  (GPC/SEC,  FTIR,  distil¬ 
lation)  that  could  not  be  settled  as  on-line  analysis  techniques,  and 
therefore  are  not  suitable  for  on-line  monitoring.  However,  use  of 
internal  standards  in  the  experiments  allowed  mass  balance  clo¬ 
sure  with  an  accuracy  of  close  to  98%,  with  the  yield  of  waxes  being 
determined  by  difference. 

A  tool  has  been  developed  in  Scilab  for  the  multivariate  analysis 
of  the  experimental  results.  The  program  developed  uses  the  pea 
function  to  determine  the  correlation  circle  and  the  scree  plot, 
and  the  fminsearch  optimization  function  for  the  orthogonal 
rotation  (varimax  rotation),  with  the  aim  being  to  facilitate  the 
interpretation  of  the  results. 

3.  Results  and  discussion 


3.1.  Experimental  results 


The  conversion  is  usually  defined  in  waste  pyrolysis  as  follows 

[39,40]: 


W0-  W 

w0-w - 


where  W,  W0  and  W x  are  the  masses  at  time  t,  at  the  beginning  and 
at  the  end  of  the  reaction.  Given  that  under  the  experimental  con¬ 
ditions  for  pyrolysis  in  the  spouted  bed  reactor,  e.g.,  short  residence 
time  and  high  N2  dilution,  no  solid  product  is  formed  in  polyolefin 
and  polyaromatic  pyrolysis,  thereby  Eq.  (1)  simplifies  to: 


W0-  W 
W0 


According  to  this  definition,  the  yield  of  each  product  or  lump 
(X,-)  is  defined  as: 


Y  -W> 

'  ■  W0 
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where  W,  is  the  mass  of  product  or  lump  i  at  a  given  time.  The  yields 
of  the  individual  products  and  lumps  obtained  in  the  thermal  pyro¬ 
lysis  of  HDPE  in  the  temperature  range  studied  in  this  paper  have 
been  published  in  a  previous  one  [41  . 

The  latter  is  the  definition  of  yield  used  in  the  multivariate 
study,  but  it  is  not  suitable  for  building  typical  X,  vs.  T  graphs. 
The  standard  scale  masks  trends  and  similar  performances  of  the 
lumps,  because  the  final  yield  that  each  lump  reaches  is  very  differ¬ 
ent.  Accordingly,  a  new  parameter  called  modified  yield,  X,*,  has 
been  defined  and  determined  as  follows: 


where  X,°°  and  W,°°  are  the  yield  and  mass  of  the  component  or 
lump  i  at  the  end  of  the  reaction,  respectively.  Fig.  1  shows  the 
results  obtained  in  the  thermal  pyrolysis  carried  out  in  batch  mode 
at  550  °C,  600  °C,  and  650  °C.  Products  have  been  grouped  according 
to  their  molecular  size  in  the  following:  CH4,  C2,  C3,  C4,  C5,  C6,  C7,  C8, 
C9,  Cio,  Cn-C2o  and  aromatics.  The  aromatics  have  been  split 
because  the  original  structure  of  the  polymer  does  not  have  any 
aromatic  rings,  so  they  are  unlikely  to  have  been  formed  as  a 
consequence  of  a  main  chain  sigma  bond  cleavage  reaction,  but 
rather  as  a  result  of  condensation  secondary  reactions  [29],  and 
they  are  coke  precursors  at  high  temperatures  [42].  A  continuous 
line  in  Fig.  1  shows  the  average  evolution  of  the  volatile  products 


(a) 


(b) 


(C) 


Fig.  1.  Evolution  with  time  of  the  modified  yield  of  the  different  lumps  grouped  according  to  their  molecular  size  in  the  batch  thermal  pyrolysis,  (a)  550  °C,  (b)  600  °C, 
(c)  650  °C. 
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(VP),  corresponding  to  the  sum  of  all  the  product  families  quantified 
by  gas  chromatography. 

At  550  °C  (Fig.  la)  methane  and  the  aromatics  have  a  similar 
trend,  with  a  slightly  slower  evolution  than  the  non-aromatic 
C5-Cn  hydrocarbons  and  the  rest  of  the  C4_  gas  components, 
which  account  for  26  wt.%  of  all  the  volatile  products  (VP)  formed. 
The  C4_  gas  components  follow  a  similar  trend  to  the  non-aromatic 
C5-Cn  hydrocarbons,  which  account  for  almost  70  wt.%  of  the  VP 
at  this  temperature.  A  wider  dispersion  is  observed  at  600  °C 
(Fig.  lb),  although  gas  components,  28  wt.%  of  the  VP,  again  follow 
a  similar  trend  (including  methane  in  the  gas  components).  How- 
ever,  remarkable  differences  are  observed  in  the  components  that 
constitute  the  gasoline  lump  (60  wt.%  VP).  C5  and  C6  components 
undergo  a  similar  formation  rate  to  the  gas  components  and  most 
of  the  VP,  but  C7  and  C8  hydrocarbons  evolve  more  slowly,  and  hea¬ 
vier  products  (C10  and  Cn-C2o)  are  formed  much  faster,  reaching 
their  maximum  yield  at  relatively  short  reaction  times.  The  exper¬ 
iments  carried  out  at  the  highest  temperature  studied,  650  °C 
(Fig.  lc),  clearly  show  the  effect  of  temperature  on  the  reaction 
kinetics,  as  the  time  needed  for  the  original  polymer  to  fully 
degrade  becomes  half  of  that  required  at  the  lowest  temperature 
studied.  A  similar  formation  rate  of  the  gas  components  (C4_)  is 
observed  again,  as  was  the  case  at  the  lower  temperatures  studied. 
These  gas  components  constitute  almost  40  wt.%  of  the  VP  at  this 
temperature.  The  formation  rate  of  these  components  is  slightly 
slower  than  those  of  liquid  components.  Among  the  liquid  compo¬ 
nents,  those  in  the  C5-C8  range  record  a  similar  trend,  whereas  the 
heaviest  ones,  C9  and  Ci0,  are  formed  faster.  These  results  are  the 
basis  for  proposing  models  in  which  waxes  are  formed  as  primary 
products  from  the  direct  cracking  of  long  polymer  chains,  and  they 
subsequently  undergo  gas  phase  secondary  cracking  reactions  to 
form  other  fractions  [29,30  . 

Fig.  2  shows  the  results  of  the  same  experiments  (HDPE  thermal 
pyrolysis  at  550  °C,  600  °C,  and  650  °C),  although  in  this  case  the 
product  components  have  been  grouped  according  to  the  nature 
of  their  C— C  bonds.  Thus,  products  have  been  grouped  as  follows: 
light  alkanes  (LA),  light  olefins  (LO),  paraffins  (P),  isoparaffins  (iP), 
olefins  (O),  benzene  (B),  toluene  (To),  ethylbenzene  +  xylenes  (Eth- 
ylBz  +  Xyl),  C9  aromatics  (C9  Arom),  Ci0  aromatics  (Ci0  Arom),  Cu 
aromatics  (Cu  Arom),  and  naphthenes.  As  in  Fig.  1,  the  average 
evolution  of  the  VP  is  shown  as  reference. 

Fig.  2a  data  (corresponding  to  550  °C)  show  that  light  olefins 
follow  a  similar  trend  to  the  VP  as  a  whole.  Light  alkanes  and  tol¬ 
uene  show  a  slower  formation  rate,  as  observed  for  methane  in 
Fig.  la.  Nevertheless,  all  the  other  aromatics  detected  at  this  tem¬ 
perature  (benzene,  ethylbenzene,  and  xylenes)  and  the  other  prod¬ 
uct  families  (paraffins,  isoparaffins,  olefins,  and  naphthenes)  have  a 
similar  formation  rate  throughout  the  reaction.  Taking  into 
account  the  properties  of  the  polymer  in  the  feed  (a  polyolefin), 
the  prevailing  products  are  components  with  double  bonds.  There¬ 
fore,  light  olefins  (C4_)  and  olefins  (C5-Cu)  account  for  more  than 
50  wt.%  of  the  VP.  The  delay  for  toluene  (the  most  common  aro¬ 
matic  at  this  temperature,  3.2  wt.%,  is  much  higher  than  the  sum 
of  the  remaining  aromatic  components)  and  light  alkanes  confirms 
that  aromatic  rings  are  formed  as  a  consequence  of  condensation 
secondary  reactions  [43  . 

Fig.  2b  corresponds  to  the  results  of  the  experiments  carried  out 
at  600  °C.  Similarly  to  the  experiments  at  550  °C,  toluene  and  the 
remaining  aromatic  components  have  a  slower  formation  rate  at 
600  °C.  The  aromatic  components  account  for  less  than  10  wt.%  of 
the  VP  at  this  temperature,  with  olefins  being  the  main  compo¬ 
nents,  40  wt.%  of  the  VP.  The  yield  of  isoparaffins  increases  in  a 
more  pronounced  way  than  at  550  °C,  reaching  35  wt.%  of  the  VP. 
Although  the  results  obtained  at  650  °C  are  not  as  disperse  as  those 
obtained  at  600  °C,  they  confirm  the  trends  observed  at  lower 
temperatures. 


3.2.  PCA  results 

Figs.  1  and  2  give  clues  about  the  common  trends  of  some  of  the 
components  and  validate  the  hypothesis  of  this  work  concerning 
the  significance  of  having  criteria  to  group  products  for  the  pro¬ 
posal  of  kinetic  schemes  and  even  for  classifying  the  products  as 
primary  or  secondary  in  order  to  progress  in  the  proposal  of  mech¬ 
anisms  (cleavage,  condensation)  for  transforming  some  products 
into  others.  Nevertheless,  the  information  obtained  from  those  fig¬ 
ures  is  somewhat  qualitative,  and  therefore  lacking  rigor.  More¬ 
over,  information  is  lost  in  these  types  of  representations  due  to 
the  dense  cloud  of  experimental  points  that  prevents  detecting 
the  time  when  maximum  production  is  achieved  (corresponding 
to  the  maximum  in  the  dX*/dt  curve).  Multivariate  analysis  tech¬ 
niques  like  PCA  are  then  planned  for  a  systematic  analysis  of  the 
trends  of  each  lump.  The  fact  that  all  the  products,  whether  they 
are  primary  or  secondary,  are  formed  from  the  raw  polymer  has 
to  be  taken  into  account  when  this  statistical  technique,  PCA,  is 
applied  to  the  study  of  pyrolysis,  because  all  the  products  are 
related  to  each  other,  as  observed  in  the  correlation  matrix. 

Another  peculiarity  of  the  application  of  PCA  to  the  batch  pyro¬ 
lysis  process,  which  has  already  been  observed  in  a  previous  paper, 
is  that  a  single  component  accounts  for  more  than  95%  of  the  total 
variance  (  fable  1).  Accordingly,  the  points  are  vertically  aligned  in 
the  correlation  circle  before  the  varimax  rotation  is  applied. 

The  results  of  the  eigenvalue  matrix  ( Table  1 )  confirm  that  the 
employment  of  only  two  components  suffices  to  explain  the  total 
variance.  Thus,  the  correlation  circle  can  be  elaborated  and  the 
varimax  rotation  applied  with  the  corresponding  2D  representa¬ 
tion  for  the  analysis  of  the  experimental  results.  Fig.  3  shows  the 
loadings  of  the  12  lumps  (variables),  in  which  the  products  of  the 
thermal  pyrolysis  of  HDPE  at  550  °C,  600  °C,  and  650  °C  have  been 
grouped  according  to  their  molecular  size  (carbon  atom  numbers). 

The  results  shown  in  Fig.  3a,  corresponding  to  the  batch  thermal 
pyrolysis  of  HDPE  at  550  °C  arranged  by  their  molecular  size,  con¬ 
firm  and  complement  the  trends  observed  in  Fig.  la.  Four  groups  or 
lumps  may  be  distinguished.  The  first  one,  as  pointed  out  in  the  X* 
vs.  t  plot,  consists  of  CH4  and  the  aromatics,  and  is  related  to  com¬ 
ponent  2.  A  second  group  is  made  up  of  C2,  C5,  and  C7  components; 
a  third  one  is  made  up  of  C3,  C4,  and  C6,  (they  stand  over  the  hor¬ 
izontal  axis  in  Fig.  3a  and  are  not  related  to  component  2);  and  a 
fourth  one  made  up  of  only  C8.  As  observed,  C8  products  are  clearly 
related  to  component  1  and  their  behavior  is  very  different  to  those 
products  that  are  solely  related  to  component  2,  i.e.,  aromatics  and 
methane.  The  ability  of  PCA  to  find  analogous  performances  is 
noteworthy.  Thus,  although  in  Fig.  la  the  performance  of  aromatics 
and  methane  is  visibly  different  to  the  remaining  products,  the 
multivariate  analysis  allows  discerning  slight  differences  among 
the  gases,  as  C2  (mainly  ethylene)  is  separated  from  C3  and  C4 
(mostly  propylene  and  butenes),  and  even  detecting  certain  differ¬ 
ences  in  the  liquid  fraction,  where  C5  and  C7  products,  mostly  ole¬ 
fins,  have  a  slightly  different  trend  compared  to  C6  and  C8  products, 
with  a  higher  concentration  of  isoparaffins.  This  information  is 
very  helpful  for  the  proposal  of  new  kinetic  schemes. 

The  correlation  circle  after  the  varimax  rotation  for  the  results 
carried  out  at  600  °C  in  batch  mode  (Fig.  3b)  shows  that,  as  in 
Fig.  lb,  the  performance  of  Ci-C6  components  is  quite  similar, 
but  different  to  that  of  C7-C9  products.  C10  and  C11+  are  the  sole 
products  that  appear  in  the  second  quadrant,  and  significantly 
affected  by  component  2.  These  results  are  explained  by  the  much 
faster  formation  of  Ci0  and  especially  Cn+  products. 

The  results  obtained  at  650  °C  when  arranging  the  products 
according  to  their  molecular  size  (Fig.  3c)  confirm  some  of  the 
groups  already  proposed  and  establish  new  ones.  All  these  results 
allow  proposing  the  lumps  that  should  be  considered  in  the  kinetic 
scheme  and  establishing  the  sequence  in  which  the  reactions  take 
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(a) 


(b) 


(c) 


Fig.  2.  Evolution  with  time  of  the  modified  yield  of  the  different  lumps  grouped  according  to  the  nature  of  their  C— C  bond  in  batch  thermal  pyrolysis,  (a)  550  °C,  (b)  600  °C, 
(c)  650  °C. 


place.  Gases  have  a  visibly  different  performance,  as  they  are  all  in 
the  third  quadrant  and  negatively  related  to  component  2,  whereas 
the  light  liquid  fraction  (C5-C8)  follows  a  different  trend.  The  per¬ 
formance  of  aromatic  components,  which  varies  with  temperature, 
suggests  they  are  secondary  products  that  are  favoured  by  an 
increase  in  temperature.  Other  components  are  also  affected  by 
secondary  reactions,  because  they  are  grouped  in  a  different  way 
when  temperature  is  increased. 

Westerhout  et  al.  [29]  proposed  that  waxes  were  the  only  pri¬ 
mary  product,  but  this  statement  has  been  subsequently  ques¬ 
tioned  by  other  authors  [30,42],  who  also  proposed  gases  and 
liquid  as  primary  products.  The  results  obtained  at  550  °C  and 


600  °C  show  that  Ci0_  products  have  a  different  performance  to 
the  heavier  hydrocarbons,  i.e.  the  latter  are  formed  sooner.  Thus, 
as  Gonzalez  et  al.  44]  proposed,  products  heavier  than  the  gasoline 
fraction,  Cu+,  should  be  considered  as  waxes  (or  pyrolysis  solid 
products)  in  the  thermal  pyrolysis  of  polyolefins  due  to  their  sim¬ 
ilar  formation  rate.  As  temperature  is  increased,  the  cracking  of  the 
polymer  is  more  severe  and  at  650  °C  even  the  lighter  fraction  of 
waxes,  CH-C20,  undergoes  secondary  cracking  reactions  to  give 
light  hydrocarbons,  Ci0_  [41].  Furthermore,  the  location  of  meth¬ 
ane  together  with  the  aromatic  components  at  550  °C,  close  to 
the  gas  products  at  600  °C,  and  the  different  trend  the  gases  follow 
at  650  °C  are  indications  that  the  gas  lump  comes  from  the  severe 
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Table  1 

PCA  initial  statistics  matrix. 


Component 

Initial  eigenvalues 

Total 

Variance  (%) 

Cumulative  (%) 

550  °C,  molecular  size 

1 

8.987 

99.85 

99.85 

2 

600  °C,  molecular  size 

0.012 

0.13 

99.98 

1 

11.746 

97.88 

97.88 

2 

650  °C,  molecular  size 

0.204 

1.70 

99.58 

1 

11.844 

98.70 

98.70 

2 

550  °C,  bond  type 

0.146 

1.22 

99.92 

1 

8.939 

99.32 

99.32 

2 

0.054 

0.60 

99.92 

600  °C,  bond  type 

1 

11.555 

96.29 

96.29 

2 

650  °C,  bond  type 

0.359 

2.99 

99.28 

1 

11.926 

99.38 

99.38 

2 

0.057 

0.47 

99.86 

cracking  of  the  different  products  that  make  up  the  wax  fraction. 
These  cracking  reactions  yield  intermediate  products  (olefins)  that 
follow  cyclization,  aromatization  reactions,  producing  aromatics 
and  light  alkanes. 

These  results  evidence  that  the  kinetic  scheme  for  HDPE  ther¬ 
mal  pyrolysis  at  the  550-650  °C  temperature  range  should  take 
into  account  at  least  5  lumps:  (i)  the  primary  products  derived 
from  the  raw  polymer  cracking,  which  are  generally  referred  to 
as  waxes;  (ii)  gaseous  products,  Ci-C4,  which  are  primary  or  sec¬ 
ondary  products  depending  on  the  reaction;  (iii)  non-aromatic 
C5-C9  hydrocarbons,  which  are  primary  or  secondary  products; 
(iv)  aromatic  products  (clearly  secondary  products);  (v)  the  char, 
which  is  a  tertiary  product  derived  from  polyaromatic  products 
[45],  although  it  has  not  been  detected  in  this  study  because  we 
operated  at  temperatures  lower  than  700  °C  and  with  relatively 
low  residence  times  for  the  pyrolysis  products.  As  a  conclusion, 
Fig.  4  shows  the  kinetic  scheme  for  HDPE  thermal  pyrolysis.  The 
continuous  arrows  are  primary  reactions,  whereas  the  dashed  ones 
are  secondary  and  tertiary  ones. 

This  scheme  is  similar  to  the  one  proposed  by  [29],  with  the 
inclusion  of  a  new  lump  (the  non-aromatic  C5-C9  liquid  fraction), 
given  that  the  formation  rate  of  these  components  cannot  be 
explained  by  any  of  the  other  lumps.  Besides,  the  scheme  does 
not  consider  the  formation  of  aromatics  from  gases  because  there 
is  no  evidence  of  this  formation  in  the  multivariate  analysis,  i.e., 
they  are  not  related  by  component  1  or  2  in  the  PCA,  and  so  an 
increase  in  one  lump  does  not  imply  a  decrease  in  the  other. 

In  order  to  decide  whether  the  criterion  of  grouping  according 
to  molecular  size  or  the  nature  of  the  C— C  bonds  is  the  more  suit¬ 
able,  Fig.  5  shows  the  results  of  the  same  experiments,  but  group¬ 
ing  the  products  into  the  following  categories:  light  alkanes  (LA), 
light  olefins  (LO),  C5-Cn  paraffins  (P),  C5-Cn  isoparaffins  (iP), 
C5-C11  olefins  (O),  benzene  (B),  toluene  (To),  and  ethylben¬ 
zene  +  xylenes  (EthylBz  +  Xyl),  C9  aromatics  (C9  Arom),  C10  aromat¬ 
ics  (C10  Arom),  Cu  aromatics  (Cn  Arom),  and  naphthenes. 

The  results  corresponding  to  the  lowest  temperature  studied 
(550  °C,  Fig.  5a)  evidence  a  similar  performance  for  most  of  the  main 
products  in  thermal  HDPE  pyrolysis;  that  is,  aliphatic  hydrocarbons, 
except  light  alkanes,  have  a  similar  formation  rate.  Light  alkanes,  as 
well  as  toluene,  are  located  in  the  second  quadrant,  being  positively 
affected  by  component  2.  This  trend  has  already  been  observed  in 
Fig.  3a  for  aromatics  and  methane.  Nevertheless,  single-ring 
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Fig.  3.  Loadings  of  the  12  lumps  (variables)  grouped  according  to  the  molecular  size 
in  the  batch  thermal  pyrolysis  of  HDPE,  (a)  550  °C,  (b)  600  °C,  (c)  650  °C. 
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Fig.  4.  Kinetic  scheme  proposed  for  the  thermal  pyrolysis  of  HDPE. 


aromatics  and  light  alkanes  are  minor  products,  as  they  account 
for  less  than  5  wt.%  of  the  product  yield  at  the  temperature  studied. 
Furthermore,  the  rest  of  the  aromatics  identified  at  this  temperature 
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differences  are  observed  for  the  aromatic  components.  Thus,  the 
lighter  ones  (benzene,  toluene,  xylenes,  ethylbenzene,  and  C9  aro¬ 
matics)  appear  grouped  in  the  graph,  but  C10  and  Cn  are  substan¬ 
tially  separated  from  the  remaining  pyrolysis  products. 

Given  that  the  yield  of  single-ring  aromatics  accounts  for  less 
than  3  wt.%  of  the  product  stream,  the  yield  of  acyclic  aliphatic 
components,  particularly  isoparaffins  and  olefins  [41  ,  is  over 
95  wt.%  and  the  information  obtained  from  PCA  by  grouping  the 
product  components  according  to  the  nature  of  the  C— C  bond 
has  not  provided  any  relevant  information  concerning  different 
performances  between  isoparaffins  and  olefins.  Therefore,  kinetic 
scheme  proposals  should  be  based  on  the  lumps  defined  according 
to  molecular  size  [23,29,30,43]  and  not  on  the  nature  of  the  C— C 
bond  [28  . 
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4.  Conclusions 

PCA  has  proven  to  be  a  versatile  method  for  verifying  the  crite¬ 
ria  used  to  group  the  large  number  of  product  components  that  are 
formed  in  the  thermal  pyrolysis  of  polyolefins.  In  fact,  PCA  allows 
the  detection  and  recognition  of  similar  performance  patterns 
among  different  compounds  or  lumps,  and  therefore  helps  to 
establish  a  precise  definition  of  lumps  and  the  hierarchy  of  the 
reactions,  differentiating  between  primary  and  secondary 
products. 

The  grouping  of  the  HDPE  thermal  pyrolysis  products,  in  the 
550-650  °C  range,  according  to  their  molecular  size,  allows  a  pre¬ 
cise  definition  of  the  components  within  the  most  common  lumps 
used  in  the  literature:  gas  (Ci-C4),  liquids  (non-aromatic  C5-Cu), 
aromatics,  and  waxes  (Cn+).  Moreover,  multivariate  analysis 
shows  that  waxes  are  the  intermediate  compounds,  and  the 
remaining  lumps  are  produced  in  secondary  reactions  from  these 
waxes  and  the  condensation  of  derived  products. 

The  grouping  of  products  according  to  bond  type  shows  that 
condensation,  cyclization  and  aromatization  producing  single-ring 
aromatics  are  clearly  affected  by  temperature.  Nevertheless,  multi¬ 
variate  analysis  applied  to  major  pyrolysis  products  (acyclic  ali¬ 
phatic  products,  especially  olefins  and  isoparaffins)  does  not 
allow  extracting  conclusive  information  for  the  proposal  of  kinetic 
schemes  based  on  the  grouping  of  compounds  according  to  hydro¬ 
carbon  families. 
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Fig.  5.  Loadings  of  the  12  lumps  (variables)  grouped  according  to  the  nature  of  their 
C— C  bond  in  the  batch  thermal  pyrolysis  of  HDPE,  (a)  550  °C,  (b)  600  °C,  (c)  650  °C. 
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(benzene,  ethylbenzene,  and  xylenes)  record  a  different  formation 
rate  compared  to  that  of  toluene.  The  difference  between  toluene 
and  the  rest  of  the  aromatics  remains  at  600  °C  (Fig.  5b),  although 
an  increase  in  temperature  induces  notable  alterations  in  the  group¬ 
ing  into  lumps.  Light  alkanes  record  a  similar  formation  rate  to  all 
the  other  aliphatic  acyclic  hydrocarbons  or  to  the  lightest  aromatic 
components  (benzene,  ethylbenzene,  and  xylenes).  Furthermore, 
the  remainder  of  the  aromatic  components  and  acyclic  hydrocar¬ 
bons  (C9— Cn  aromatics,  and  naphthenes)  are  located  in  the  fourth 
quadrant,  being  closely  related  to  component  1. 

Finally,  the  correlation  circle  of  the  results  obtained  at  650  °C 
when  the  components  are  arranged  according  to  the  nature  of  their 
C— C  bonds  (Fig.  5c)  corroborates  that  the  formation  rate  of  all  the 
acyclic  aliphatic  hydrocarbons  is  similar.  Moreover,  the  main 


References 

[1]  Plastics  Europe,  Plastics  -  the  Facts,  2013;  An  analysis  of  European  latest 
plastics  production,  demand  and  waste  data. 

[2]  G.  Dodbiba,  N.  Haruki,  A.  Shibayama,  T.  Miyazaki,  T.  Fujita,  Combination  of 
sink-float  separation  and  flotation  technique  for  purification  of  shredded  PET- 
bottle  from  PE  or  PP  flakes,  Int.  J.  Miner.  Process.  65  (2002)  11-29. 

[3]  P.J.  Donaj,  W.  Kaminsky,  F.  Buzeto,  W.  Yang,  Pyrolysis  of  polyolefins  for 
increasing  the  yield  of  monomers’  recovery,  Waste  Manage.  32  (2012)  840- 
846. 

[4]  S.M.  Al-Salem,  S.  Evangelisti,  P.  Lettieri,  Life  cycle  assessment  of  alternative 
technologies  for  municipal  solid  waste  and  plastic  solid  waste  management  in 
the  Greater  London  area,  Chem.  Eng.  J.  244  (2014)  391-402. 

[5]  R.  Aguado,  M.  Olazar,  B.  Gaisan,  R.  Prieto,  J.  Bilbao,  Kinetics  of  polystyrene 
pyrolysis  in  a  conical  spouted  bed  reactor,  Chem.  Eng.  J.  92  (2003)  91-99. 

[6]  S.H.  Jung,  S.J.  Kim,  J.S.  Kim,  Fast  pyrolysis  of  a  waste  fraction  of  high  impact 
polystyrene  (HIPS)  containing  brominated  flame  retardants  in  a  fluidized  bed 


364 


R.  Aguado  et  al./ Chemical  Engineering  Journal  254  (2014)  357-364 


reactor:  The  effects  of  various  Ca-based  additives  (CaO,  Ca(OH)(2)  and  oyster 
shells)  on  the  removal  of  bromine,  Fuel  95  (2012)  514-520. 

[7]  A.  Adrados,  I.  de  Marco,  B.M.  Caballero,  A.  Lopez,  M.F.  Laresgoiti,  A.  Torres, 
Pyrolysis  of  plastic  packaging  waste:  a  comparison  of  plastic  residuals  from 
material  recovery  facilities  with  simulated  plastic  waste,  Waste  Manage.  32 
(2012)  826-832. 

[8]  G.  Elordi,  M.  Olazar,  G.  Lopez,  M.  Amutio,  M.  Artetxe,  R.  Aguado,  J.  Bilbao, 
Catalytic  pyrolysis  of  HDPE  in  continuous  mode  over  zeolite  catalysts  in  a 
conical  spouted  bed  reactor,  J.  Anal.  Appl.  Pyrolysis  85  (2009)  345-351. 

[9]  A.  Lopez-Urionabarrenechea,  I.  de  Marco,  B.M.  Caballero,  M.F.  Laresgoiti,  A. 
Adrados,  A.  Aranzabal,  Catalytic  pyrolysis  of  plastic  wastes  with  two  different 
types  of  catalysts:  ZSM-5  zeolite  and  Red  Mud,  Appl.  Catal.  B  104  (2011)  211  — 
219. 

[10]  A.  Coelho,  L.  Costa,  M.M.  Marques,  I.M.  Fonseca,  M.  Lemos,  F.  Lemos,  The  effect 
of  ZSM-5  zeolite  acidity  on  the  catalytic  degradation  of  high-density 
polyethylene  using  simultaneous  DSC/TG  analysis,  Appl.  Catal.  A  413  (2012) 
183-191. 

[11]  G.  Elordi,  M.  Olazar,  P.  Castano,  M.  Artetxe,  J.  Bilbao,  Polyethylene  cracking  on  a 
spent  FCC  catalyst  in  a  conical  spouted  bed,  Ind.  Eng.  Chem.  Res.  51  (2012) 
14008-14017. 

[12]  H.  Bockhorn,  A.  Hornung,  U.  Hornung,  D.  Schwaller,  Kinetic  study  on  the 
thermal  degradation  of  polypropylene  and  polyethylene,  J.  Anal.  Appl. 
Pyrolysis  48  (1999)  93-109. 

[13]  A.  Aboulkas,  K.  El  Harfi,  A.  El  Bouadili,  Thermal  degradation  behaviors  of 
polyethylene  and  polypropylene.  Part  I:  pyrolysis  kinetics  and  mechanisms, 
Energy  Convers.  Manage.  51  (2010)  1363-1369. 

[14]  G.  Elordi,  M.  Olazar,  G.  Lopez,  M.  Artetxe,  J.  Bilbao,  Continuous  polyolefin 
cracking  on  an  HZSM-5  zeolite  catalyst  in  a  conical  spouted  bed  reactor,  Ind. 
Eng.  Chem.  Res.  50  (2011)  6061-6070. 

[15]  S.M.  Al-Salem,  P.  Lettieri,  J.  Baeyens,  Recycling  and  recovery  routes  of  plastic 
solid  waste  (PSW):  a  review,  Waste  Manage.  29  (2009)  2625-2643. 

[16]  J.  Ceamanos,  J.F.  Mastral,  A.  Millera,  M.E.  Aldea,  Kinetics  of  pyrolysis  of  high 
density  polyethylene.  Comparison  of  isothermal  and  dynamic  experiments,  J. 
Anal.  Appl.  Pyrolysis  65  (2002)  93-110. 

[17]  A.  Durmus,  S.N.  Koc,  G.S.  Pozan,  A.  Kasgoz,  Thermal-catalytic  degradation 
kinetics  of  polypropylene  over  BEA,  ZSM-5  and  MOR  zeolites,  Appl.  Catal.  B  61 
(2005) 316-322. 

[18]  J.M.  Encinar,  J.F.  Gonzalez,  Pyrolysis  of  synthetic  polymers  and  plastic  wastes. 
Kinetic  study,  Fuel  Process.  Technol.  89  (2008)  678-686. 

[19]  S.M.  Al-Salem,  P.  Lettieri,  Kinetic  study  of  high  density  polyethylene  (HDPE) 
pyrolysis,  Chem.  Eng.  Res.  Des.  88  (2010)  1599-1606. 

[20]  A.  Coelho,  L.  Costa,  M.D.  Marques,  I.  Fonseca,  M.A.  Lemos,  F.  Lemos,  Using 
simultaneous  DSC/TG  to  analyze  the  kinetics  of  polyethylene  degradation- 
catalytic  cracking  using  HY  and  HZSM-5  zeolites,  React.  Kinet.  Mech.  Catal.  99 
(2010)  5-15. 

[21]  B.  Saha,  P.K.  Reddy,  A.K.  Ghoshal,  Hybrid  genetic  algorithm  to  find  the  best 
model  and  the  globally  optimized  overall  kinetics  parameters  for  thermal 
decomposition  of  plastics,  Chem.  Eng.  J.  138  (2008)  20-29. 

[22]  S.S.  Kim,  S.  Kim,  Pyrolysis  characteristics  of  polystyrene  and  polypropylene  in  a 
stirred  batch  reactor,  Chem.  Eng.  J.  98  (2004)  53-60. 

[23]  F.  Ding,  L.  Xiong,  C.R.  Luo,  H.R.  Zhang,  X.D.  Chen,  Kinetic  study  of  low- 
temperature  conversion  of  plastic  mixtures  to  value  added  products,  J.  Anal. 
Appl.  Pyrolysis  94  (2012)  83-90. 

[24]  J.  Bilbao,  A.T.  Aguayo,  J.M.  Arandes,  Coke  deposition  on  silica-alumina 
catalysts  in  dehydration  reactions,  Ind.  Eng.  Chem.  Prod.  Res.  Dev.  24  (1985) 
531-539. 

[25]  W.  Kaminsky,  Possibilities  and  limits  of  pyrolysis,  Makromol.  Chem. 
Macromol.  Symp.  57  (1992)  145-160. 


[26]  J.A.  Conesa,  A.  Marcilla,  R.  Font,  Kinetic-model  of  the  pyrolysis  of  polyethylene 
in  a  fluidized-bed  reactor,  J.  Anal.  Appl.  Pyrolysis  30  (1994)  101-120. 

[27]  R.  Aguado,  M.  Olazar,  B.  Gaisan,  R.  Prieto,  J.  Bilbao,  Kinetic  study  of  polyolefin 
pyrolysis  in  a  conical  spouted  bed  reactor,  Ind.  Eng.  Chem.  Res.  41  (2002) 
4559-4566. 

[28]  W.C.  Huang,  M.S.  Huang,  C.F.  Huang,  C.C.  Chen,  K.L.  Ou,  Thermochemical 
conversion  of  polymer  wastes  into  hydrocarbon  fuels  over  various  fluidizing 
cracking  catalysts,  Fuel  89  (2010)  2305-2316. 

[29]  R.W.J.  Westerhout,  JAM.  Kuipers,  W.P.M.  van  Swaaij,  Experimental 
determination  of  the  yield  of  pyrolysis  products  of  polyethene  and 
polypropene.  Influence  of  reaction  conditions,  Ind.  Eng.  Chem.  Res.  37  (1998) 
841-847. 

[30]  G.  Elordi,  G.  Lopez,  M.  Olazar,  R.  Aguado,  J.  Bilbao,  Product  distribution 
modelling  in  the  thermal  pyrolysis  of  high  density  polyethylene,  J.  Hazard. 
Mater.  144  (2007)  708-714. 

[31]  E.  Jarde,  F.  Vilmin,  L.  Mansuy,  P.  Faure,  PY-GC/AED  and  chemometric 
correlation  to  characterize  sewage  sludges  of  different  origins,  J.  Anal.  Appl. 
Pyrolysis  71  (2004)  553-567. 

[32]  D.  Meier,  I.  Fortmann,  J.  Odermatt,  O.  Faix,  Discrimination  of  genetically 
modified  poplar  clones  by  analytical  pyrolysis-gas  chromatography  and 
principal  component  analysis,  J.  Anal.  Appl.  Pyrolysis  74  (2005)  129-137. 

[33]  P.  Khare,  B.P.  Baruah,  P.G.  Rao,  Application  of  chemometrics  to  study  the 
kinetics  of  coal  pyrolysis:  a  novel  approach,  Fuel  90  (2011)  3299-3305. 

[34]  Z.  Sebestyen,  Z.  May,  K.  Reczey,  E.  Jakab,  The  effect  of  alkaline  pretreatment  on 
the  thermal  decomposition  of  hemp,  J.  Therm.  Anal.  Calorim.  105  (2011)  1061- 
1069. 

[35]  R.  Aguado,  A.  Arrizabalaga,  M.  Arabiourrutia,  G.  Lopez,  J.  Bilbao,  M.  Olazar, 
Principal  component  analysis  for  kinetic  scheme  proposal  in  the  thermal  and 
catalytic  pyrolysis  of  waste  tyres,  Chem.  Eng.  Sci.  106  (2014)  9-17. 

[36]  G.  Elordi,  M.  Olazar,  R.  Aguado,  G.  Lopez,  M.  Arabiourrutia,  J.  Bilbao,  Catalytic 
pyrolysis  of  high  density  polyethylene  in  a  conical  spouted  bed  reactor,  J.  Anal. 
Appl.  Pyrolysis  79  (2007)  450-455. 

[37]  G.  Elordi,  M.  Olazar,  M.  Artetxe,  P.  Castano,  J.  Bilbao,  Effect  of  the  acidity  of  the 
HZSM-5  zeolite  catalyst  on  the  cracking  of  high  density  polyethylene  in  a 
conical  spouted  bed  reactor,  Appl.  Catal.  A  415  (2012)  89-95. 

[38]  M.  Arabiourrutia,  G.  Elordi,  G.  Lopez,  E.  Borsella,  J.  Bilbao,  M.  Olazar, 
Characterization  of  the  waxes  obtained  by  the  pyrolysis  of  polyolefin 
plastics  in  a  conical  spouted  bed  reactor,  J.  Anal.  Appl.  Pyrolysis  94  (2012) 
230-237. 

[39]  G.  Lopez,  R.  Aguado,  M.  Olazar,  M.  Arabiourrutia,  J.  Bilbao,  Kinetics  of  scrap  tyre 
pyrolysis  under  vacuum  conditions,  Waste  Manage.  29  (2009)  2649-2655. 

[40]  M.  Amutio,  G.  Lopez,  R.  Aguado,  J.  Bilbao,  M.  Olazar,  Biomass  oxidative  flash 
pyrolysis:  auto  thermal  operation,  yields  and  product  properties,  Energy  Fuels 
26  (2012)  1353-1362. 

[41]  G.  Elordi,  M.  Olazar,  G.  Lopez,  M.  Artetxe,  J.  Bilbao,  Product  yields  and 
compositions  in  the  continuous  pyrolysis  of  high-density  polyethylene  in  a 
conical  spouted  bed  reactor,  Ind.  Eng.  Chem.  Res.  50  (2011)  6650-6659. 

[42]  N.  Horvat,  F.T.T.  Ng,  Tertiary  polymer  recycling:  study  of  polyethylene 
thermolysis  as  a  first  step  to  synthetic  diesel  fuel,  Fuel  78  (1999)  459-470. 

[43]  R.W.J.  Westerhout,  J.  Waanders,  JAM.  Kuipers,  W.P.M.  van  Swaaij,  Recycling  of 
polyethene  and  polypropene  in  a  novel  bench-scale  rotating  cone  reactor  by 
high-temperature  pyrolysis,  Ind.  Eng.  Chem.  Res.  37  (1998)  2293-2300. 

[44]  Y.S.  Gonzalez,  C.  Costa,  M.C.  Marquez,  P.  Ramos,  Thermal  and  catalytic 
degradation  of  polyethylene  wastes  in  the  presence  of  silica  gel,  5A  molecular 
sieve  and  activated  carbon,  J.  Hazard.  Mater.  187  (2011)  101-112. 

[45]  M.  Guisnet,  P.  Magnoux,  Organic  chemistry  of  coke  formation,  Appl.  Catal.  A 
212  (2001)  83-96. 


